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Amendments to the Claims; 

This listing . of claims will replace all prior versions, and 
listings, of claims in the application: 

Listing of Claims; 

1. (Original) A method performed by a computer for 
generating a geometric pattern from an image having a plurality 
of ridges and mesh points, the method comprising the steps of: 

forming a partial differential equation by transferring 
values for positions in the image to corresponding coefficients 
of the partial differential equation; 

determining simultaneous difference equations corresponding 
to the partial differential equation and the image mesh points; 

solving the simultaneous difference equations; and 

mapping the solutions of the simultaneous difference 
equations to respective positions on the image to determine 
features of the image. 



2. (Original) The method of claim 1, wherein the image is 
a fingerprint. 

3. (Original) The method of claim 1, wherein the image is 
a facial image. 

4. (Original) The method of claim 1, wherein the image is 
a hand-palm image. 

5. (Original) The method of claim 1, wherein the image is 
an eye iris image. 
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6. (Original) The method of claim 1, wherein the image is 
a texture image. 

7. (Original) The method of claim 1, wherein the image is 
an eye retina image. 

8. (Original) The method of claim 1, wherein the step of 
forming a partial differential equation comprises the steps of: 

calculating a plurality of intrinsic properties of the 

image according to image ridge patterns- 
mapping the plurality of intrinsic properties into 

coefficients of the partial differential equation; and 
determining a boundary condition for the partial 

differential equation from the image to establish a relationship 

between properties of the image and the partial differential 

equation. 

9. (Original) The method of claim 1, wherein the step of 
determining simultaneous difference equations comprises the step 
of replacing each derivative and variable items of the partial 
differential equation with respect to each mesh point in the 
image . 

10. (Original) The method of claim 1, wherein the step of 
forming a partial differential equation comprises the steps of: 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

11. (Original) The method of claim 10, wherein the step of 
determining initial conditions comprises the steps of: 
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normalizing the image to reduce variations in gray-level 

values along ridges of the image; 

estimating property values of the image; and 
mapping the estimated property values into weight 

coefficients of the partial differential equation. 



12. (Original) The method of claim 10, wherein the step of 
determining boundary condition comprises the steps of: 
drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 



13. (Original) The method of claim 11, wherein the step of 
normalizing comprises the steps of: 

determining M, mean of the gray- level in a region R by 
computing 

M=(1/N)ZZ F(I,J) 

aJ)eR (2a) 
where N is total number pixels in the region R; F(I, J) is 
gray value of the image at point (I, J); 

determining V, variance of the region R by computing 

V=( 1/N) Z I (F(IJ)-M(F))*(F(U>M(F)) 

(^'J)^^ (2b); and 



determining a normalized region R by computing 

Ra, J) = m + sqrt((v * (F(I, J) - M) * (F(I, J) - M))A^, ifa, J) > M; ^ ^ ^ ^ 
R(I, J) = m - sqrt((v * (F(I, J) - M) * (F(I, J) - M))AO, otherwise ^ 

where m and v are the desired mean and variance values, 
respectively. 
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14. (Currently Amended) The method of claim 11, wherein 
the step of estimating property values comprises the steps of: 
dividing a region R into blocks of size b*b as B(k); 
confuting gradients at each pixel in R as 



where (I, J) and (I, J) are the gradient magnitude in x 
and y directions at pixel (I, J) of the image respectively, pi, 
p2 are positive numbers, [[p2]] £3 is negative number, and d is 
a constant expressed as step of the gradients p = pl+p2+p3; 

estimating local orientation of each block (B(k)) centered 
at pixel (I, J) by computing 



ex(I, J) = (pi * F(I - d, J) + p2 * F(I, J) + p3 * (F(I + d, J))/p, 
5x(I, J) = (pi * F(I, J - d) + p2 * F(I, J) + p3 * (Fa, J + d))/p 



(3a) 



4x(IJ)=IZ2*ax(u,v)*ay^(uy), 

(u,v)GB(k) 
^ y(I,J) = IZ(aax(u,v)*5ay(u,v)), 

(u,v)6B(k) 
e(I,J) = (l/2)atan{^(I, J)/0 ^ yI,J)} 



(3b) 



(3c) 



(3d) 



Where 6(1,}) is an estimate of the local ridge 
orientation at the block centered at pixel (I, J) ; 



computing ridge oriental vector as 



P = (Vn)Zcos(2*0(I,J)) 



(4a) 
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Q = (l/«)Zsin(2*^(/,J)) 



(4b) 



where P and Q are first and second components of the 
ridge oriental vector, respectively, and n is the total number 
of pixels calculated at the local region R; 

for each block centered at pixel (I, J), computing the 
minimal value and maximal value the pixel; 

for each block centered at pixel (I, J), computing a 
sequence of pixels seq(I, J) that take minimal and maximal value 
along the direction (a, b) , where (a, b) is orthogonal vector of 
the main oriental vector (p, q) ; 

computing freq(I, J), frequency of seq(I, J) at each block 
centered at pixel (I, J) according to the differential value 
between connected elements in seq(I, J); and 

estimating a local ridge frequency W by computing 

W = (l/K)ZIfreq(u,v) 
(u,v)e wnd(I,J) 

15. (Original) The method of claim 11, wherein the step of 
mapping comprises the step of determining weight coefficients, 
Al, A2, A3, A4, A5 and A6 of a partial differential equation 




dXdY 



^.a'u ^.au ^.du . 

+ A3 — r + A4 — + A5 — + A6*U = 0 
dY^ ax dY 



(1) 



by computing 

A1 = P*P*(P*P + Q*Q)*W*W, 

A2 = 2*(sqrt(u*u-P*P*W*W)*sqrt(v-Q*Q*W*W))AV, 
A3 = Q*Q(P*P + Q*Q)*W*W, 



{6a) 
(6b) 
(6c) 
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A4 = u*q + v, 

A5 = -v*p-u, and 

A6 = a*(P*P + Q*Q) + b 



(6d) 
(6e) 
(6f) 



where a, b, u, v are constants. 

16. (Original) The method of claim 12, wherein the step of 
setting boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 



where S is a continuous boundary defined on the discrete 
boundary B, and Bl (x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the boundary S, 
respectively. 

17. (Original) The method of claim 1, wherein the step of 
determining simultaneous difference equations comprises the 
steps of: 

integral i zing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on the mesh points for 
numerating the partial differential equation and a boundary 
condition; and 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations . 

18. (Original) The method of claim 17, wherein the step of 
integralizing comprises the steps of: 



U 



= Bl(x,y), au=B2(x,y) 
(x,y)€S as (x,y)eS 



(7) 
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denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, coitputing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, I = 0, 1, 2, , W(F), (8a) 

Y(J) = YO + J*H, J = 0, 1, 2, , H(F), (8b) 

where, (XO, YO) is top left point of the image, W(F) is the 
width of the image and H{F) is the height of the image. 

19, (Currently Amended) The method of claim 17, wherein 
the step of discretizing comprises the steps of: 

computing the derivatives in the partial differential 
equation 

Al — r+A2 + A3 — r + A4 — + A5 — + A6*U = 0 (1) 

ax' axaY av' ax ax 



with respect to each inner mesh point as 

^«[U(X + H,Y)-UpC,Y)]/H (9a) 
^«[U(X,Y+H)-U(X,Y)]/H (9b) 
~[U(X i -II,Y)-2*U(X,\0 + U(X-II,\0]/II*II) 

ax 
aHj 
ax' 



[U(X + H, Y) - 2 * U(X, Y) + U(X - H, Y)]/(H * H) ( 9c ) 



^ « [U(X, Y + H) - 2 * U(X, Y) + U(X, Y - H)]/(H * H) ( 9d) 

dY 

«[U(X + H,Y + H)-U(X + H,Y)-U(X,Y + H) + U(X,Y)]/(H*H) (9e) 



axax 

where (X, Y) is inner mesh point in region R, IMP(R); 
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discretizing the boundary condition 

^=Bl(x,yX aU 
(x,y)eS 5s 

combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

20. (Original) The method of claim 19, wherein the step of 

discretizing the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 

condition with discrete function F(I, J), wherein (I, J) is 

inner mesh point of a region in the images- 
replacing the continuous function Bl (x, y) in the boundary 

condition with a numerical function according to 
D 1 (X, Y) = n * F(X, Y) + £2, (X, Y) € IB (10) 

where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, Y) 
is the gray value at point (X, Y) on the integral boundary IB; 
and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X, Y) = n * [F(X1, Yl) - F(X, Y)]/h, (X, Y) e IB 

h = sqrt((Xl-X)*(Xl-X) + (Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 



= B2(x,y) 
(x,y)eS 



and 



(7) 



-9- 



Appln No. 10/047,535 

Amdt date September 20, 2005 

Reply to Office action of August 8, 2005 

21. (Original) The method of claim 1, wherein the step of 
mapping the solutions comprises the steps of: 

for each region computing minimum value min (R) and 

maximum value max (R) in the solutions- 
computing ratio value r (R) = 255/ (max (R) - min (R) ) ; 
for each point (I, J) in the region R, assigning v(I, J) * 

3 as the solution at point (I, J) ; 

computing W(I, J) = r(R)*(v(I, J)-min{R)) as the mapping of 

w(I, J) into gray level byte at the position (I, J) ; and 

enhancing the region R is by placing the value W{I, J) at 

position (I, J) . 

22 . (Canceled) 

23- (Currently Amended) The A metho d of claim 22 performed 
by a computer for extracting features from an image, the method 
comprising the steps of; 

establishing a mathematical model according to regional 

conditions in the image; 

converting the mathematical model into numerical equations; 

solving the numerical equations; and 

mapping respective values of the solutions of the numerical 

equations to respective regions of the image , wherein the step 
of establishing a mathematical model comprises the steps of: 

forming a partial differential equation; 

calculating a plurality of intrinsic properties of the 
image according to image ridge pattern; 

mapping the plurality of intrinsic properties into 
coefficients of the partial differential equation; and 

determining a boundary condition for the partial 
differential equation from the image to establish a relationship 
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between properties of the image and the partial differential 
equation. 



24. (Canceled) 

25. (Canceled) 

26. (Currently Amended) The method of claim [[24]] 22 , 
wherein the step of determining boundary condition comprises the 
steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 



27-29. (Canceled) 



30. (Original) The method of claim 26, wherein the step of 
setting boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 



U 



= Bl(x,y), aU 
(x, y) e S 5s 



= B2(x,y) 
(x,y)€S 



(7) 



where S is a continuous boundary defined on the discfete 
boundary B, and Bl (x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the boundary S, 
respectively. 



31-34. (Canceled) 
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35. (Currently Amended) The method of claim [[22]] 23 , 
wherein the step of tranof erring mapping respective values of 
the solutions comprises the steps of: 

for each region R in the image, computing minimum value min 

(R) and maximum value max (R) in the solutions- 
computing ratio value r (R) = 255/ (max (R) - min (R) ) ; 
for each point (I, J) in the region R, assigning v(I, J) * 

3 as the solution at point (I, J) ; 

computing W(I, J) = r(R)*(v(I, J)-min(R)) as the mapping of 

w(I, J) into gray level byte at the position (I, J) ; and 

enhancing the region R is by placing the value W(I, J) at 

position (I, J) . 

36. (Currently Amended) The method of claim [[22]] 23 , 
wherein the image is one or more of a fingerprint image, a 
facial image, a hand-palm image, an eye iris, an eye retina, and 
a texture image. 

37. (Original) A digital signal processor (DSP) having 
stored thereon a set of instructions including instructions for 
generating geometric pattern from an image having a plurality of 
ridges and mesh points, when executed, the instructions cause 
the DSP to perform the steps of: 

forming a partial differential equation by transferring 
values for positions in the image to corresponding coefficients 
of the partial differential equation; 

determining simultaneous difference equations corresponding 
to the partial differential equation and the image mesh points; 

solving the simultaneous difference equations; and 

mapping the solutions of the simultaneous difference 
equations to respective positions on the image to determine 
features of the image. 
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38. (Original) The DSP of claim 37, wherein the step of 
forming a partial differential equation comprises the steps of: 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

39. (Original) The DSP of claim 38, wherein the step of 
determining initial conditions comprises the steps of: 

normalizing the image to reduce variations in gray-level 

values along ridges of the images- 
estimating property values the image; and 
mapping the estimated property values into weight 

coefficients of the partial differential equation. 

40. (Original) The DSP of claim 38, wherein the step of 
determining boundary condition comprises the steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 

41. (Original) The DSP of claim 39, wherein the step of 
normalizing comprises the steps of: 

determining M, mean of the gray- level in a region R by 
computing 

M = (l/N)EZFa,J) 
(I,J)eR 

where N is total number pixels in the region R; F(I, J) is 
gray value of the image at point (I, J) ; 

determining V, variance of the region R by computing 
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V = (1/N)ZZ(F(I.J)-M(F))*(F(I,J)-M(F)) 

(2b) ; . 

(I,J)eR 

determining a normalized region R by computing 



R(I, J) = m + sqrt((v * (F(I, J) - M) * (F(I, J) - M))A^),if (I, J) > M; 
R(I, J) = m - sqrt((v * (F(I, J) - M) * (F(I, J) - M))A^), otherwise 



(2c) 



where m and v are the desired mean and variance values, 
respectively. 

42. (Currently Amended) The DSP of claim 39, wherein the 
step of estimating property values comprises the steps of: 
dividing a region R into blocks of size b*b as B(k.); 
computing gradients at each pixel in R as 



where (I, J) and (I, J) are the gradient magnitude in x 
and y directions at pixel (I, J) of the image respectively, pi, 
p2 are positive numbers, [ [p2] ] p3 is negative number, and d is 
a constant expressed as step of the gradients p = pl+p2+p3; 

estimating local orientation of each block (B(k)) centered 
at pixel (I, J) by computing 



ax(I, J) = (p 1 * F(I - d, J) + p2 * F(I, J) + p3 * (F(I + d, J))/p, 
ay(I, J) = (p 1 * F(I, J - d) + p2 * F(I, J) + p3 * (F(I, J + d))/p 



(3a) 



^(I, J) = EE2* 5x(u,v)*ay(u,v). 
(u,v)€B(k) 

Cy(i,J) = ll(aax(u,v)*aay(u,v)), 

(u,v)GB(k) 
^(I,J) = (l/2)atan{^(I, J)/^y(I,J)} 



(3b) 



(3c) 



(3d) 
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Where is an estimate of the local ridge orientation 
at the block centered at pixel (I, J) ; 



computing ridge oriental vector as 



P = (l/n)Zcos(2*^(I,J)) 

(IJ)€R 
Q = (l/n)Isin(2*^(I,J)) 

(I,J)eR 



(4a) 



{4b) 



Where P and Q are first and second components of the 
ridge oriental vector, respectively, and n is the total number 
of pixels calculated at the local region R; 

for each block centered at pixel (I, J), computing the 
minimal value and maximal value the pixel; 

for each block centered at pixel (I, J) , computing a 
sequence of pixels seq(I, J) that take minimal and maximal value 
along the direction (a, b) , where (a, b) is orthogonal vector of 
the main oriental vector (p, q) ; 

computing freq(I, J), frequency of seq(I, J) at each block 
centered at pixel (I, J) according to the differential value 
between connected elements in seq(I, J); and 

estimating a local ridge frequency W by computing 



W = (l/K)SSfreq(u,v) 
(u,v)€ wnd(I,J) 



(5) 
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43. (Original) The DSP of claim 39, wherein the step of 
mapping comprises the step of determining weight coefficients, 
Al, A2, A3, A4, A5 and A6 of a partial differential equation 



Al — :r + A2 + A3 — ^ + A4 — + A5 — + A6*U = 0 

ax' axav av' ax av 



(1) 



by computing 

Al = P*P*(P*P-f Q*Q)*W*W, (6a) 
A2 = 2*(sqrt(u*u-P*P*W*W)*sqrt(v-Q*Q*W*W))AV, (6b) 

A3 = Q*Q(P*P + Q*Q)*W*W, (6c) 

A4 = u*q + v, (6d) 

A5 = -v*p-u, and (6e) 

A6 = a*(P*P + Q*Q) + b (6f) 
where a, b, u, v are constants. 

44. (Original) The DSP of claim 40, wherein the step of 
setting boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 



U 



= Bl(x,y), aU 
(x,y)GS as 



= B2(x,y) 
(x,y)€S 



(7) 



where S is a continuous boundary defined on the discrete 
boundary B, and Bl (x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the boundary S, 
respectively. 
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45. (Original) The DSP of claim 37, wherein the step of 
determining simultaneous difference equations comprises the 
steps of: 

integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; . 

discretizing the image based on the mesh points for 
numerating the partial differential equation and a boundary 
condition; and 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations . 

46. (Original) The DSP of claim 45, wherein the step of 
integralizing comprises the steps of: 

denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, I = 0, 1, 2, , W(F), (8a) 

Y(J) = YO + J*H, J = 0, 1, 2, , H(F), (8b) 

where, (XO, YO) is top left point of the image, W(F) is the 
width of the image and H(F) is the height of the image. 

47. (Currently Amended) The DSP of claim 45, wherein the 
step of discretizing comprises the steps of: 

computing the derivatives in the partial differential 
equation 

a'u ^^a'u ^,au ^^au ^ 

Al — ;- + A2 + A3 — :r + A4 — + A5 — + A6*U = 0 (1) 

ax' dXdY ay' ax av 

with respect to each inner mesh point as 
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— «[U(X + H,Y)-U(X,Y)]/H (9a) 
5X 

— « [U(X, Y + H) - U(X, Y)]/H ( 9b ) 
5Y 

~[U(X I H,Y)-2*U(X,Y) + U(X-II,\0]/II*n) 

5X 
5'U 



«[U(X + H,Y)-2*U(X,Y) + U(X-H,Y)]/(H*H) (9c) 



[U(X,Y + H)-2*U(X,Y) + U(X,Y-H)]/(H*H) (9d) 



ay' 
a'u 



« [U(X + H, Y + H) - U(X + H, Y) - U(X, Y + H) + U(X, Y)]/(H * H) ( 9e ) 



axax 

where (X, Y) is inner mesh point in region R, IMP(R); 
discretizing the boundary condition 



U 



= Bl(x,y), aU 
(x,y)6S as 



= B2(x,y) 

; and ( 7 ) 

(x,y)eS 



combining the numerical derivatives of the partial 
differential equation and numerical boundary condition , 

48. (Original) The DSP of claim 47, wherein the step of 
discretizing the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F(I, J), wherein (I, J) is 
inner mesh point of a region in the image; 

replacing the continuous function Bl (x, y) in the boundary 
condition with a numerical function according to 

D1(X, Y) = fl* F(X, Y) + f2, (X, Y) € IB (10) 
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where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, Y) 
is the gray value at point (X, Y) on the integral boundary IB; 
and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X, Y) = fl*[F(Xl,Yl) .F(X,Y)]/h, (X, Y) € IB ^ ^ 

h = sqrt((Xl-X)*(Xl-X) + (Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

49. (Original) The DSP of claim 37, wherein the step of 
mapping the solutions comprises the steps of: 

for each region R, computing minimum value min (R) and 
maximum value max (R) in the solution; 

computing ratio value r (R) = 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J) ; 

computing W(I, J) = r(R)*(v(I, J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J); and 

enhancing the region R is by placing the value W{I, J) at 
position (I, J) . 

50. (Original) The DSP of claim 37, wherein the image is 
one or more of a finger print, a facial image, a hand-palm 
image, an eye iris, an eye retina, and a texture image. 

51-55. (Cancelled) 
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56. (Currently Amended) The A method of claim 51 performed 
by a computer for biometric image processing, the DSP comprising 
the steps of: 

establishing numerical relationship between visual 
appearance of the biometric image and a partial differential 
equation model of the image; and 

approximating solutions of the partial differential 

equation with a boundary condition according to the established 
numerical relationship to determine features of the biometric 
image , wherein the step of approximating solutions comprises the 
steps of: 

integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on mesh points for numerating 
a partial differential equation and a boundary condition; 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations; 

replacing the continuous function U{X, Y) in the boundary 
condition with discrete function F(I, J), wherein (I, J) is 
inner mesh point of a region in the image; 

replacing the continuous function Bl (x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y) = fl*F(X,Y) + f2, (X,Y)eIB (10) 

where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, Y) 
is the gray value at point (X, Y) on the integral boundary IB; 
and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 
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D2(X,Y) = fl*[F(Xl,Yl)-F(X,Y)]/h, (X,Y)€lB ^^^^ 

h = sqrt((Xl-X)*(Xl-X) + (Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

57, (Original) The method of claim 56, wherein the step of 
integralizing comprises the steps of: 

denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, I = 0, 1, 2, , W(F), (8a) 

Y(J) = YO + J*H, J = 0, 1, 2, , H(F), (8b) 

where, (XO, YO) is top left point of the image, W(F) is the 
width of the image and H{F) is the height of the image. 

58. (Currently Amended) The method of claim 56, wherein 
the step of discretizing comprises the steps of: 

computing the derivatives in the partial differential 
equation 

Al — r+A2 + A3 — 7 + — + A5 — + A6*U = 0 (1) 

aX' dXdY dY' dX BY 

with respect to each inner mesh point as 

^ « [U(X + H, Y) - U(X, Y)]/H ( 9a) 

— « [U(X, Y + H) - U(X, Y)]/H ( 9b ) 

9Y 
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~[U(X I II,\0-2*U(X,\0 + U(X-II.'^0]/II*II) 

dX 



8X' 



« [U(X + H, Y) - 2 * U(X, Y) + U(X - H, Y)]/(H * H) (9c) 



^ « [U(X, Y + H) - 2 * U(X, Y) + U(X, Y - H)]/(H * H) ( 9d) 

dY 

« [U(X + H, Y + H) - U(X + H, Y) - U(X, Y + H) + U(X, Y)]/(H * H) { 9e ) 

dXdX 

where (X, Y) is inner mesh point in region R, IMP (R) ; 



discretizing the boundary condition 



U 



= Bl(x,y), aU 



(x,y)€S 3s 



= B2(x,y) 

; and (7) 

(x,y)€S 



combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

59. (Original) The method of claim 56, wherein the step of 

discretizing the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 

condition with discrete function F(I, J), wherein (I, J) is 

inner mesh point of a region in the images- 
replacing the continuous function Bl (x, y) in the boundary 

condition with a numerical function according to 
D1(X, Y) = fl * F(X, Y) + G, (X, Y) € IB (10) 

where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, Y) 
is the gray value at point (X, Y) on the integral boundary IB; 
and 
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replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X,Y) = fl*[F(Xl,Yl)-F(X,Y)]/h, (X,Y)€lB ^^^^ 

h = sqrt((Xl-X)*(Xl-X) + (Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

60. (Original) The method of claim 56, further comprising 
the step of mapping the solutions of the partial differential 
equation to respective positions on the image to determine 
features of the image. 

61. (Original) The method of claim 60, wherein the step of 
mapping the solutions comprises the steps of: 

for each region R, computing minimum value min (R) and 
maximum value max (R) in the solution; 

computing ratio value r (R) = 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J) ; 

computing W(I, J) = r(R)*(v(I, J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J); and 

enhancing the region R is by placing the value W(I, J) at 
position (I, J) . 

62. (Currently Amended) The method of claim [[51]] 56, 
wherein the image is one or more of a finger print, a facial 
image, a hand-palm image, an eye iris, an eye retina, and a 
texture image. 
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